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Recent work on single bubble sonoluminescence (SBSL) has shown that many features of 
this phenomenon, especially the dependence of SBSL intensity and stability on experi- 
mental parameters, can be explained within a hydrodynamic approach. More specifically, 
many important properties can already be derived from an analysis of bubble wall dy- 
namics. This dynamics is conveniently described by the Rayleigh-Plesset (RP) equation. 
In this work we derive analytical approximations for RP dynamics and subsequent ana- 
lytical laws for parameter dependences. These results include (i) an expression for the 
onset threshold of SL, (ii) an analytical explanation of the transition from diffusively 
unstable to stable equilibria for the bubble ambient radius (unstable and stable sono- 
luminescence), and (iii) a detailed understanding of the resonance structure of the RP 
equation. It is found that the threshold for SL emission is shifted to larger bubble radii 
and larger driving pressures if surface tension is enlarged, whereas even a considerable 
change in liquid viscosity leaves this threshold virtually unaltered. As an enhanced vis- 
cosity stabilizes the bubbles against surface oscillations, we conclude that the ideal liquid 
for violently collapsing, surface stable SL bubbles should have small surface tension and 
large viscosity, although too large viscosity {rji > AOrjwater) will again preclude collapses. 



1. Introduction 

1.1. Sonoluminescence 

The analysis of the dynamics of a small bubble or cavity in a fluid dates back to the 
work of Lord Rayleigh (1917) at the beginning of this century. A large number of 
publications followed in subsequent decades, including the studies of oscillating bubbles 
by Plesset (1949, 1954), [Eller fc Crum (1970)| , Flynn (1975a, 1975b), [Lauterborn (1976) 



Prosperetti (1977), Plesset & Prosperetti (1977), and others. In recent years, a renascence 



of bubble dynamics has occurred initiated by the discovery of single bubble sonolumi- 



nescence (SBSL) by Gaitan (1990), see also Gaitan et al. (1992) 
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SBSL is an intriguing phenomenon: A single gas bubble of only a few /xm size, levi- 
tated in water by an acoustic standing wave, emits light pulses so intense as to be visible 
to the naked eye. The standing ultrasound wave of the driving keeps the bubble in 
position at a pressure antinode and, at the same time, drives its oscillations. The ex- 
periments of Putterman's group (Barber & Putterman 1991; Barber et al. (1994,1995); 
Hiller et al. 1994; Lofstedt, Barber & Putterman 1993; Lofstedt et al. 1995; Weninger, 
Putterman & Barber 1996) and others have revealed a multitude of interesting facts 
about SBSL: the width of the light pulse is small (Barber & Putterman 1991 give 50 ps 
as upper threshold, Moran et al. 1995 10 - recent measurements by Gompf et al. 1997 
report 100-300ps, depending on the forcing pressure and gas concentration in the liq- 
uid), the spectrum shows no features such as lines (Hiller, Putterman & Barber 1992; 
Matula et al. 1995). While the exact mechanism of light emission is still an open is- 
sue, almost all suggested theories - see e.g. Lofstedt et al. (1993), Hiller et al. (1992), 



Fhnt fc Sushck (1989)i |Wu fc Roberts (1993)| , |Frommhold fc Atchley (1994)| , Moss et al. (1994), 
Bernstein & Zakin (1995), Moss, Clarke & Young (1997) - agree that temperatures of 
at least 10'*-10^ K are reached during bubble collapse. This, together with the light in- 
tensity, clearly shows that SBSL relies on an extraordinarily powerful energy focusing 
process. 



In our previous publications 


Brenner, Lohse & Dupont (1995), Brenner et al. (1996a, 


1996b), Hilgenfeldt, Lohsc & Brenner (1996) 




Brenner, Hilgenfeldt & Lohse (1996), 


Lohse et al. (1997) 


, and Lohse & Hilgenfeldt (1997) 


we calculated phase diagrams for 



bubbles and have focused on the identification of parameter regimes where SBSL occurs. 
As a scan of the whole multi-dimensional parameter space is by far too expensive for full 
numerical simulations of the underlying fundamental equations (i.e., Navier-Stokes and 
advection-diffusion PDEs), it is necessary to introduce approximations. The necessary 
conditions for SL to occur could be calculated from the dynamics R{t) of the bubble wall, 
which is - apart from a tiny interval around the bubble collapse - very well described by 
the Rayleigh-Plesset (RP) equation. We call this approach the RP-SL bubble approach. 

The key parameters in an SL experiment are the ambient bubble radius Rq (ra- 
dius under normal conditions of 1.013 x 10^ Pa = latm and 20°C), the driving pres- 
sure amplitude Pa, and the gas concentration in the water surrounding the bubble 
Pod/Pq, measured by its partial pressure divided by the ambient pressure. Note that 
Rq is not at the experimenter's disposal, but adjusts itself by gas diffusion on a slow 
time scale of seconds. Its size can, however, be measured in experiment, e.g., by Mie 



scattering techniques as in Barber et al. (1995) or by direct microscopic imaging, see 



Tian, Ketterling fc Apfel (1996) , Holt fc Gaitan (1996) . On time scales much smaller 
than those of diffusive processes, e.g. for one period of driving, Rq may be regarded as a 
constant to high accuracy. 



In Hilgenfeldt et al. (1996) we found that the Pa/ Po^Poo/ Po state space is divided into 
regions where (diffusively) stable SL, unstable SL or no SL are to be expected, in excellent 
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agreement with experimental findings. These results will now be briefly presented in the 
following subsection. 



1.2. Stability requirements 

Stable sonoluminescence is characterized by light emission in each period of driving at 
precisely the same oscillation phase and precisely the same brightness for millions (and 
sometimes billions) of cycles. We found that it occurs in a tiny section of the whole 
parameter space only, and that the calculated domain agrees very well with experimen- 
tal findings, cf. Hilgenfeldt et al. (1996), Lohsc et al. (1997). Its boundaries are set by 
certain dynamical and stability conditions imposed upon the oscillating bubble (Brenner 
et al. 1995, Brenner et al. 1996a, Hilgenfeldt et al. 1996): (i) The bubble wall velocity 
during collapse must reach the speed of sound in the gas Cg to ensure sufficient energy 
transfer from the liquid to the gas. (ii) The bubble must be stable towards non-spherical 
oscillations of its surface which lead to fragmentation. Bubble fragments have meanwhile 
been experimentally observed by J. Holzfuss (private communication, 1997). (iii) The 
bubble must be stable towards diffusive processes, i.e., it must not dissolve or grow by 
rectified diffusion; diffusively growing bubbles show unstable SL. A further requirement 
of (iv) chemical stability becomes important when the bubble contains molecular gases 
which are able to dissociate and recombine with liquid molecules (Brenner et al. 1996, 
Lohse et al. 1997). E.g., the differences in the parameter regimes of SL in air bubbles vs. 
SL in noble gas bubbles can consistently be accounted for by dissociation of N2 and O2 
in an air bubble; these molecular constituents of air are burned, leaving only inert gases 
in the bubble (the experimental work of Holt & Gaitan 1996 supports this model). We 
therefore restrict ourselves for simplicity to the case of a bubble filled with argon. An 
extension to reactive gas mixtures as analysed in Lohse et al. (1997) is straightforward. 
Also, we specify the liquid in which the bubble oscillates to be water, as in most SBSL 
experiments. 

Figure |l| illustrates how the conditions (i)-(iii) determine domain boundaries in the 
Pa~Ro state space. Criterion (i) means that the Mach number with respect to Cg is larger 
than 1, i.e., 

|M,| = ^>1, (1.1) 

and it is fulfilled for bubbles with large enough ambient radii Rq and at large enough 
forcing amplitude Pa, i.e., right of the dashed line in figure The shape stability 



condition (ii) - see picsset (1954)| , [Birkhoff (1954)| , |Eller fc Crum (1970)| , [Strube (1971 



Prosperetti (1977), Brenner et al. (1995) or Hilgenfeldt et al. (1996) for detailed stud- 



ies -, on the other hand, limits the parameter domain in which bubbles can stably 
oscillate to small i?o ^ 4 — 5 fim, within our boundary layer approximation. As the RP 
SL approach neglects effects of thermal conduction, which has a damping influence on 
surface oscillations, this upper limit on Rq may be somewhat higher in reality. Holt & 
Gaitan's (1996) experimental results seem to give a threshold around 7 fim. (i) and (n) 
together determine the shaded area of potentially sonoluminescing bubbles in figure 0. 
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Figure 1. Stability conditions for a bubble in Pa-Ro parameter space. Bubbles above the 



line \Mg\ = 1 fulfill the energy focusing condition (1.1). Bubbles below the shape instability 
lines are stable towards non-spherical surface oscillations. The solid line represents the (long 
time scale) parametric instability, the dashed line short time scale shape instabilities, for de- 
tails see Hilgenfeldt et al. (1996). Bubbles on the thick lines are in diffusive equilibrium for 
Poo/Po ~ 0.2,0.002, respectively. Thick arrows indicate regions of bubble growth and shrinking 
by diffusive processes. 

The actual position of a stable SL bubble in Pa-Ro parameter space is determined by 
condition (iii) for stable diffusive equilibria of the gas inside the bubble and the dissolved 
gas in the liquid (thick lines in figure ^) . These equilibrium lines show negative slope 
whenever the equilibria are unstable, i.e., bubbles below the line shrink and dissolve, 
bubbles above the line grow. At large gas concentrations in the liquid (e.g. Poo/Pq ^ 0.2, 
left curve), only unstable equilibria are possible in the parameter range of interest. Tiny 
ratios Poo /^b ~ 0.002 (right curve) are necessary for diffusive stability (i.e., the fluid must 
be strongly degassed) . The positive slope of the upper branch of the curve characterizes 



these bubbles as stable. The computation of diffusive equilibria is explained in § 2.2 



1.3. Summary of results of the present work 

Having identified the parameter regions for SBSL through numerically solving the RP 
equation, the question arises if one can understand the shape and size of these regions 
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analytically, i.e., by analysing the bubble dynamics equations. In principle, all of the 
conditions that determine the occurrence of stable/unstable/no SL depend only on prop- 
erties of bubble dynamics. Therefore, we set out in this work to derive analytical approx- 
imations for RP dynamics and subsequently find scaling laws or approximate analytical 
expressions for our numerical curves presented above, in order to give a clearer insight 
into the role of different physical processes governing the dynamical equations. More- 
over, more practical reasons make analytical expressions highly desirable, as the multi- 
dimensional parameter space of SBSL experiments cannot be scanned in detail just by 
numerical solution of the RP equation. Our analytical efforts strongly build on previous 
work, most notably that of Lofstedt et al. (1993). We present the most important results 
in this subsection, written such that experimental parameters can be directly inserted to 
yield numerical values. Here we have used fixed w = 27r x 26.5 kHz and Pq = 1 atm. More 
detailed results and the complete derivations for general uj, Pq will be given in the corre- 
sponding Sections. All the presented approximations naturally have limited parameter 
regimes of validity, which include the region of sonoluminescing bubbles in all cases. 

We will demonstrate in §H that, in order to understand the location of diffusive equi- 
libria, it is sufficient to analyse the parameter dependence of the ratio of the maximum 
bubble radius to its ambient radius {Rmax/Ro), see Lofstedt et al. (1993). We show in 
§ ^ that two clearly distinct kinds of bubble dynamical behaviour exist depending on 
and Rq: weakly oscillating and strongly collapsing bubbles. The transition between these 
two states is rather abrupt and occurs for given Pa at an ambient radius 

" 9^ Pa-Po Pa/ Pa -I ^ ' 

This transition is controlled by the surface tension a, i.e., strong collapses are easier to 
achieve for small cr. 

In we derive analytical approximations to RP dynamics for all phases of the os- 
cillation cycle of a strongly collapsing bubble. We find that in this regime the bubble 
essentially collapses like an empty cavity (see Rayleigh (1917)) according to 

p \ 3/5 / 1^ (\ 2/5 



with the time of maximum bubble compression t* and the driving period T = 2i:/lo. 

Following the collapse, a series of characteristic afterbounces of the bubble radius 
occurs. We show in § 4.3 that they are the cause for the wiggly structure of the diffusive 
equilibrium curves and the \Mg\ = 1 line in figure |l|. The location of the wiggles can 
be understood as a parametric resonance phenomenon. A Mathieu approximation yields 
the ambient radius of the k^^ wiggle as 

R\, « 37.0 Mm • ^J^jf=^l + 0.487 Mm • [-^JT^ + ) (1-4) 

with the abbreviation g = (1 + Pal Pa)- 

Section 4.4 deals with the bubble expansion. In the regime of strong bubble collapses, 
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an approximate result for the dependence of the maximum radius on Pa and Rq is 

« 67.2 + 0.112 ( ^) ' + 99.5(PjFo - . (1.5) 
fim \fj,m J 

With Rmaxi the location of diffusive equilibria in {Pa,Ro) parameter space can be cal- 
culated. 

A closer discussion of the role of surface tension and viscosity of the liquid rji is pre- 
sented in §1^. In particular, the viscosity of water is so small that it has no significant 
influence on bubble dynamics. Oscillations are only viscosity-dominated if 

/ 0.487 um\ „„^fRo\ , 

which corresponds to r]i > 'iOrjwater for typical Rq. 

Note that these equations are not fit formulas, but are analytically derived from the 
RP dynamics. They are all verified by comparison to full numerical solutions in the 
appropriate domains of validity. With these formulas, we are able to understand most of 
the parameter dependences of SL analytically. Section ^ presents conclusions. 



2. Rayleigh Plesset bubble dynamics 

2.1. Notation and parameters 

Since Lord Rayleigh (1917, see Lamb 1932 for earlier references) treated the collapse of 
an empty cavity in a liquid, a lot of refinement has been done in the modelling of the 
dynamics of spherical domain walls in liquids. The main step towards bubble dynamics 
was the introduction of a variable external driving pressure and of the influence of surface 
tension by Plesset (1949) . 

An ODE for the bubble radius can be derived from the Navier-Stokes equations from 
an approximation valid to the order of R/ci, where R is the speed of the bubble wall 
and Q is the sound speed in the liquid. Following Prosperetti & Lczzi (1986), Lofstedt 
et al. (1993) and many others, we will henceforth denote the following ODE as Rayleigh- 
Plesset (RP) equation: 



pi (^RR + = Pgas{R, t) - P{t) - Po 



^ d , , i? 2cr 

+ --Pgas{R,t)-Ar]i-~—. 2.1 
ci at R R 

The left-hand side of this ODE for the bubble radius R consists of dynamical pressure 

terms already known to Rayleigh {pi — 1000 kgm~^ is the density of water). Pq = 1 atm 

is the constant ambient pressure, P{t) the ultrasound driving, modelled as a spatially 

homogeneous, standing sound wave, i.e., 

P{t) = ~Pa COS Ujt = -Po ■ p cos LUt (2.2) 

with the dimensionless forcing pressure amplitude p = Pa/ Pq and a fixed frequency of 
tt" = 27r X 26.5 kHz (period T « 38 ps), which is a common value in many experiments 
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like those of Barber et al. (1994) and Hiller et al. (1992). The wavelength of this sound 
in water is about 5 cm, while the bubble radii treated in this work never exceed 200 /iw. 
Because of this separation of scales, it is common to assume spatial homogeneity, as 
stated above. We will refer to the sum of experimentally controllable pressures as the 
external pressure p^xt = -Po + Pit)- By definition, the external pressure exerts maximally 
outward directed forces {pext — Poi^ — p) < ^) on the bubble at t = 0. 

The other terms on the right-hand side of equation (2.1) model the influence of the 
surface tension at the bubble-water interface (cr = 0.073 kg s^^), the water viscosity 
(rji — 1.00 X 10^"^ Pas), and of emitted sound waves from the bubble (cf. Keller & Miksis 
1980, this term contains the speed of sound in water q ~ 1481 to s^^). 

The gas pressure pgas{R, t) inside the bubble is assumed to obey a van der Waals type 
process equation 



PgasiR,t) =PgasiR{t)) 



2a 

Rq 



(2.3) 



Rq being the ambient bubble radius and h the (collective) van der Waals hard core ra- 
dius h = i?o/8.86 (for argon) (Lide 1991). The pressure exerted by surface tension was 
included explicity in (2.1). The a dependence of the prefactor of the polytropic expres- 
sion ensures that Rq is the radius of a static (unforced) bubble, neglecting effects of gas 
diffusion. Note that (2J) presupposes homogeneity of the pressure inside the bubble. 
This is of course not satisfied in the final stages of bubble collapse, as a more detailed 



investigation of the gas dynamics inside the bubble reveals, cf. Wu fc Roberts (1993' 



Moss et al. (1994)| , [Vuong fc Szeri (1996)| , [Evans (1996")| , |B renncr et 



(1996b)| , Moss 

et al. (1997), but the violent collapse phase lasts only ^ Ins out of the T « 38 /zs 
of the oscillation cycle. Therefore, this approximation does not severely affect our analy- 
sis of bubble wall dynamics. We furthermore set the effective polytropic exponent k w 1 
as for this frequency and bubble ambient radii below 20 /im the bubbles can be consid- 
ered to be isothermally coupled to the surrounding liquid (Plesset & Prosperetti 1977), 
except during the small time interval around the bubble collapse, where the extremely 
rapid bubble dynamics requires adiabatic treatment of the gas. This will be taken into 
account in § 4J and 4^. The solid line of figure ^(a) shows a time series R{t) from (2.1) 
for relatively strong driving Pa = l.Aatm and moderate ambient radius i?o — 4.0 /ito. 
The typical feature of the oscillations of R{t) is a slow expansion for approximately half 
a cycle of driving, followed by a rapid and violent collapse and a series of afterbounces 
corresponding to an almost free oscillation of the bubble. The time scale of the after- 
bounces is thus set by the period of the bubble's (small amplitude) eigenoscillations, 
whose frequency We ~ 1 MHz can be easily obtained from a linearization of (^.l[ ) : 

9 , . Pn cos LUt 

R + uil{R-Ro)^-^-— — with (2.4) 



piRo 



3Po 



where we have set k = 1 and neglected surface tension and viscosity effects. Including 
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Figure 2. (a) Bubble dynamics for Pa — lAatm, Ro = 4.0 /im resulting from the RP equation 
(2.1) (solid), Flynn's equation (Al) (dotted) and Gilmore's equation (A2) (dashed). The inset 
shows a blowup of the vicinity of the collapse, (fe) External pressure pext for the dynamics in 
(a), tm is the time of maximum expansion, t* the time of collapse; Pext — a,t t = t+,t-. The 
different intervals of the oscillation cycle treated in Section ^ are indicated, (c) Mach numbers 
Ml for the time interval displayed in the inset of (a). For the RP equation (solid) and Flynn's 
equation (dotted) Mi = R/ci with constant c;, for Gilmore's equation (dashed) Mi = R/Ci with 
pressure dependent C;. 



surface tension yields 

, SPn 4(T / 2 > , , , 
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pressure term 


definition 


Pace 


piRR 


Pvel 




Pgas 




Psur 


R 


Pvis 




Psnd 




Pext 


Po — Pa COS Ljt = Po{l — p COSut) 



Table 1. Definition of tfie pressure terms in tfie RP equation (2.1) used in tliis work. 



where as = 2a / {PqRq) is the ratio of surface tension pressure to Pq at R — Rq. as ~ 1 
for Rq w l.S/iTO, while for larger Rq it becomes very small. 

The resonance radius, on the other hand, is defined as the ambient radius of a bubble 
with uje = ^, i-e., 

/ op \ 1/2 

Rres = ^ ~ lObfim . (2.7) 



For convenience, we list in table |^ the definition of the different pressure terms of ( ^.l| 
which will appear throughout this paper. 



Besides the solution of the RP equation (2.1), figure U{a) shows time series obtained 



from other commonly used bubble dynamical equations, namely Flynn's and Gilmore's 
equation, which are discussed in detail in Appendix A. It is obvious that, for bubbles 
in the SBSL regime, all equations yield very similar R{t) dynamics. It is only upon 
magnification of the small time interval around the collapse (figure ^) that the differences 
between these descriptions of bubble dynamics becomes apparent. 

The deviations of the RP, Flynn, and Gilmore equations from each other may be- 
come pronounced when the bubble is driven at very high pressure amplitudes such as 
Pa — 5 aim (cf. Lastman & Wentzell 1981). These pressures are common in cavitation 
fields, but they are far too high to allow for stable bubbles in SBSL experiments (with 
the possible exception of SBSL in high magnetic fields described by Young, Schmiedel 
& Kang 1996). Sonoluminescent bubbles require a driving pressure amplitude in a nar- 
row window 1.1 atm ^ Pa 1.5 atm. It is this range of Pa that we will mainly focus 



on in this work. Only in § 4.4 results in the range of cavitation field pressures will 



briefly be displayed. Direct and indirect measurements of the size of SL bubbles e.g. in 



Barber fc Putterman (1992)| , Tian et al. (1996), or |Hoh fc Gaitan (1996)| indicate that 



typical i?o he around 5 /im. 

2.2. Calculating diffusive eguilibria from RP dynamics 

A computation of points of diffusive equilibrium in the Pa-Ro plane from first princi- 
ples requires solution of an advection diffusion PDE with appropriate boundary condi- 
tions, coupled to the RP equation. This is numerically far too expensive to allow for 
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a scan of the whole Pa-Ro parameter space. In 



Brenner et al. (1996a) and Hilgenfeldt 



et al. (1996), we therefore employed the method introduced by ^yrillas fc Szeri (1994) 



and Lofstedt et al. (1995) , which is based on the separation of the driving time scale T 
and the diffusive time scale Tdiff S> T. Within this approximation, the task is massively 
reduced to the solution of the RP equation and the computation of weighted averages of 
the form 

(/). = £ fmHt)dt j £ RHt)dt . (2.8) 

The mass flux into or out of the bubble is then proportional to Poo — (j>gas) ^ (see Fyrillas 
& Szeri 1994). An equilibrium point is characterized by the simple condition 



Pea — {Pgas)^ 



and it is stable if 



d {Pga. 



s/4 



dRo 



(2.9) 



(2.10) 



is positive. 

Figure displays {pgas)/^ for different Pa. The graphs show characteristic wiggles 



for larger Rq (which can be explained from resonance effects, see §4.3) and, for large 
enough Pa, a global minimum at some critical Rq — Rq. If Rq > Rq, even with no 
wiggles present, the bubbles are diffusively stable according to the sign of the slope 
(3. For small Rq < Rq, all equilibria are unstable, i.e., the bubble either dissolves or 



grows by rectified diffusion, see Blake (1949), Filer & Flynn (1964) (the latter case can 



lead to unstable SBSL, cf. Hilgenfeldt et al. 1996). The possibility of multiple stable 



equilibria because of the resonance structure was recognized earlier by Church (1988) and 



Kamath, Prosperetti fc Fgolfopoulos (1993)| . Here we analyse the formal and physical 



origin of the positive overall slope of {pgas) ^ (Ro) for large Rq, which is an essential 
property of stable SBSL bubbles. 

In the average {pgas) ^ the pressure is weighted with i?^(t) and will therefore be domi- 
nated by the value oipgas at Rmax- For large radii, we can neglect the excluded volume 



in the van der Waals formula and replace (2.2) by an ideal gas law under isothermal 
conditions. 



{Pga 



gas/ 4^ 



Pn 



(l + a^) 



RlR{t)At 



Rq 



Rn 



(2.11) 



^ is a prefactor that is due to the different shape of the integrands R{t) and R'^{t). 
A crude estimate of £, can be obtained by approximating R(t) by a parabola R(t) ~ 
Rmaxi^ - IGt^/T"^) and integrating R and R"^ over one half cycle from -T/4 to T/A. This 
gives ^ — 105/64 sa 1.64, which is quite accurate in reproducing numerical results. With 
this saddle point approximation introduced by Lofstedt et al. (1993), the key parameter 
for diffusive equilibria is the expansion ratio Rmax/Ro- Figure |^ demonstrates the close 
relation between {pgas) 4^ and Rmax/Ro as functions of Rq. The expansion ratio displays 
a maximum at Rq, corresponding to the minimum of {pgas)^^. In order to determine 
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Figure 3. (a) Average pressure (pgas)^- The curves are parametrized with the dimensionless 
driving pressures p = 1.0, 1.1, 1.2, 1.3, 1.4, and 1.5 from top to bottom. Note the logarithmic 
scale of the ordinate. The horizontal dashed line is an example for a given gas concentration 
Poo/Po in the liquid. Crosses mark the corresponding _Ro where diffusive equilibria are found, 
(b) Expansion ratio Rmax/Ro for the same p as in (a): p — 1.0 to 1.5 in steps of 0.1 from bottom 
to top. The maxima occur at almost exactly the same _Ro ~ Rq as the minima of {pgas)^- Up 
arrows mark the Rq, down arrows indicate the transition radii Rq^ for p=1.3, 1.4, 1.5. 



diffusive equilibrium points, one has to look for the intersections of the {pgas)^ / Po curves 
in figure||with a horizontal line given by poo/-Po (cf. equation (2.9)). Note that degassing 
to tiny partial pressures is necessary to achieve equilibria in the Rq range of pure argon 



SL bubbles; this fact was first realized by Lofstedt et al. (1995) 



For high enough Pa, there are two equilibrium values for i?o, the larger one being 
a stable equilibrium, the smaller one being unstable. If Pa is decreased, {Pgas) ^ / Pq 
increases and the equilibria come closer together. This can also be seen in figure |: for 
decreasing Pq, the Rq values given by the Pao/Po = 0.002 equilibrium curve approach 
each other. Eventually, at a certain Pa the stable und the unstable equilibrium coalesce 
and for smaller Pa no equilibrium is possible. This is reflected in figure || by the fact that 
the whole {pgas) ^ / Pq curve lies above Poo/Po- 

For relatively high gas concentrations such as Pac/Po = 0.2, stable equilibria can only 
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exist for very large Rq, where the bubbles are shape unstable. But if the concentration is 
lowered, e.g. to Poo/Po — 0.002, the stable branch (positive slope in figure |l|) enters the 
region of sonoluminescent bubbles, whereupon stable SL can set in. The occurrence of 
stable and unstable branches depends on the existence of a minimum in {pgas)^, which 
in turn necessitates a maximum in Rmax/ Rq (figure and b). Therefore, to analyse 
the lines of diffusive equilibria in figure 0, it is sufficient to explain the maximum of 
the expansion ratio figure 0(5) and its dependence on _Ro and p; this question will be 



addressed in S 4.4 



3. Quasistatic Blake threshold 

The transition from sharply increasing Rmax/ Rq for small Rq to decreasing expansion 
ratios for large Rq (figure ^) marks an important boundary between two very different 
types of bubble dynamics. Consider figure ^ where two examples of bubble dynamics 
for the same Pa — 1.5 atm and only minutely different ambient radii are displayed. The 
smaller bubble exhibits a weak (although obviously not sinusoidal) oscillation with a 
maximum expansion ratio Rmax/Ro ~ 2; no collapse is visible. The time series of the 
larger bubble is almost indistinguishable from the other until i « 0. But then, a rapid 
expansion to Rmax/Ra ~ 10 occurs, followed by a strong collapse, the typical dynamics 
of a sonoluminescing bubble, cf. figure §(a). 

Figure H shows the compression ratio Rmin/Ro of the minimum radius achieved during 
bubble oscillation to the ambient radius as a function of Pa and Rq ■ A sharp transition, 
like in the expansion ratio, is obvious in this graph and it occurs at the same Rq. For small 
Pa and small Rq, Rmin/Ro is near one; we denote such bubbles as weakly oscillating. For 
large Pa and Rq, a horizontal plane at Rmin/Ro ~ /i/i?o indicates collapse to a radius 
very near the hard core radius. We say that these latter bubbles exhibit strong collapses. 

The key to understanding this transition from weakly oscillating to strongly collapsing 
bubbles lies in the existence of a threshold for spontaneous bubble expansion known as 
the Blake threshold (Blake 1949, Atchley 1989). It is normally considered for bubbles 
under static conditions: let us first set Pa (and thus also Pext) constant in time, and 
correspondingly take R{t) to be time- independent. Then the RP equation reduces to 

where for pgas again the isothermal ideal gas law was used, which is certainly an excellent 
approximation for the static situation. For p^xt > 0, equation (3.1) has exactly one 
solution for positive R, and it corresponds to a stable equilibrium. If pext < but small 
in absolute magnitude, two equilibria exist, the one at larger R being unstable, i.e., a 
bubble with larger radius would grow indefinitely. Finally, at a critical pf^^. < (Blake 
threshold pressure, cf. Prosperetti 1984) the two equilibrium points merge and disappear 
in an inverse tangent bifurcation. In this situation, pgas is always larger than pe^t +Psur 
and (|3.l[) cannot be fulfilled for any radius. Thus, the assumption of a time-independent 
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Figure 4. Bubble dynamics for Pa = 1.5 aim, Ro = 0.85 /im (dashed) and Rq = 0.90 /ijn 
(solid). The larger bubble undergoes dynamical expansion and a strong collapse, the smaller 
oscillates weakly with the driving field (lower part). 



R{t) has to be dropped. A dynamical expansion ensues with significant contributions 
from the dynamical pressure terms on the left-hand side of (^?l|). 

Returning to the oscillatory driving pi^xt = ^o(l — pcosojt), we notice that the driving 
period T — 2tt / oj k, jis is long compared to the time scale of the bubble's eigenoscilla- 
tions 27r/we ~ 1 /is. Thus, we can consider the external pressure oscillations as quasistatic 
and follow Blake's argument as above. As Pext < is necessary to cross the Blake thresh- 
old, we must require p > 1 here. Obviously, the most sensitive point in the cycle is t = 0, 
where Pext is negative and of magnitude (p — l)^o- 

The quasistatic approximation (|3.1[) describes the complete time series of a weakly 



oscillating bubble with good accuracy. Rewriting (3.1), we obtain the cubic equation 

{pcosLot-l)R''-'^R^+(l + ^\Rl^O. (3.2) 
^0 V naff) J 

Given a time t for which Pf-^t < 0, there is a critical Rq — R^ above which the two 
positive real solutions of (3.2) become complex. When this happens, the weak oscillation 
dynamics is no longer a valid description and the transition to strong collapses occurs. 
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Figure 5. Compression ratio Rmin/Ro as a function of Pa and -Ro- The two regimes of bubble 
dynamics are clearly visible: weakly oscillating bubbles for small -Ro and small Pa, strong 
collapses to the hard core radius for large -Ro and large Pa- 



For given p, the smallest transition radius Rq' is required for t ~ 0. For R^J, therefore, 



the discriminant of (3.2) at t — must vanish, i.e. 



3 2cr 2 32 cr^ 
^° + 7^^°~27P3b-l)2 ^^■^> 

After a lengthy but straightforward calculation, the transition ambient radius R^ at 
given p = Pa I Pa is 

— X 1/3 



Rtr^l^ ( 2 2 ^ 

3Po]\{p-ir ^ip-i)yip-ir 




(3.4) 



Note that R^ is a real number for all p. In figure ^ the calculated Rq" from (3.4) is 
compared to the numerical values (identified by the condition Rmin/Ro = 0.5). The 
agreement is very good, the errors at higher Pa being only about 0.01 /im. 

When Rq exceeds -Rg'', there is a period of time around t — where the right-hand 



side of (3.1) cannot be zero, but must be positive. Then, the dynamical terms neglected 
so far must become noticeable and a dynamical expansion follows which can only be 
stopped when p^xt has again become large enough to allow for a stable radius equilibrium. 
When the bubble growth is stopped, the expanded bubble does not experience significant 
outward directed forces and, consequently, undergoes a violent collapse. If Rq is only 
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slightly larger than Rq^ , the time scale separation still holds for a large portion of the 
cycle, cf. figure |[ 
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Figure 6. Transition ambient radii Rq^ from numerical solution of the RP ODE (circles) and 
from (3.4) (dashed). This figure shows the same parameter range as figure ^, from which the 
\Mg\ = 1 curve was taken (thin dashed line). The transition to collapsing bubbles occurs for 
slightly smaller pressures and ambient radii than the onset of light emission at \AIg \ — 1. 



It is immediately obvious from ( |3.3[ ) and that surface tension plays a key role 

in this transition mechanism from weak oscillations to strong collapses. If p > 1, weak 
oscillations at small Rq are dominated by the influence of cr, whereas strongly collapsing 
(larger) bubbles are controlled by the properties of dynamical expansion and collapse 
(cf. §i 4T|, 4.4). Note that in a fluid with very small cr, already bubbles of very small size 
will show collapses (see also Lofstedt et al. 1995 and Akhatov et al. 1997). It should also 
be emphasized here that the crucial driving parameter for the transition is {p — 1), i.e., 
the difference of driving pressure amplitude Pa and ambient pressure Pq, rather than Pa 
itself. 



In the limit of large forcing p ^ 1, (3.3) yields the much simpler formula 



Tjtr 



Pop- 



(3.5) 
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It can be seen from figure ^ that in this limit the difference between Rq^ (onset of 
transition) and R^ (extremum of expansion and compression ratio) becomes neghgibly 
small. Thus, (3^) is also a good approximation to the critical Rq we were trying to 
identify. This is confirmed by figure ^ from which also the (small) errors of the saddle 
point approximation (determining Rq from R,nax/Ro instead of (^305)4) can be read off. 




Figure 7. Critical ambient radii Rq from numerical computation of the extrema of ipgaa)^ 
(filled circles) and Rmax/Ro (open diamonds) and from the asymptotic law for R!^ equation 



( ]3.5| ) (dashed line). The scaling behaviour Rq ex. l/{p — 1) from (3.5) is quite accurate for higher 
pressures p > 1.3. 



What is the maximum radius of a bubble weakly oscillating at R\^(j))l Inserting (3.5) 



into (3.2) with t — Q and expanding to the same order in l/(p — 1) gives 



Rmax = (3.6) 

3 Po P - 1 

in the large p limit. This yields a minimum expansion ratio of Rmax/Ro' ~ for the 
onset of bubble collapse, which is an analytical justification of Flynn's (1975b) definition 
of a transient cavity. In that work, a strongly collapsing bubble was characterized by an 
expansion ratio > 2. 

As the collapse sets in rather abruptly when i?o is enlarged, we expect that R^q also 
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marks the transition to bubbles which fulfill the Mach criterion (1.1). Figure || shows the 
\Mg\ — 1 line of figure |l] together with the R([{p) line according to (3^). Both curves 
display the same trend, approaching each other at large p. The Blake transition occurs 
for smaller Pa and Rq than those necessary for \Mg\ > 1, i.e., for possible light emission. 
The physical consequence of this is that, upon increasing the driving force, the bubble 
first emits cavitation noise due to collapses and only afterwards starts to emit light. 
Indeed, such a sequence of events has been reported by W. Eisenmenger & B. Gompf 
(private communication, 1996). 

The transition line Rq^{p) is shifted towards smaller Rq for smaller a. This means that 
collapses of the same violence can be achieved (for a given Rq range) with smaller driving 
pressures in a liquid with less surface tension. Note however that such bubbles will also 
be stronger affected by surface instabilities, whereas in a liquid with high a, bubbles are 
more surface stable. It is therefore possible to obtain violent collapses at larger Rq in 
liquids with larger surface tension using larger driving pressures. 



4. A guided tour of RP dynamics 

Let us now explain in detail the dynamics of strongly collapsing bubbles (as shown e.g. 
in figure ||a) . To this end, we divide the oscillation cycle of the bubble into several time 
intervals indicated in figure |^(6), where is the time of maximum bubble radius, t* the 
time of minimum bubble radius (after collapse), and t+ = —t^ = arccos(l/p)/w the time 
when pext changes its sign from positive (contracting) to negative (expanding) values. 
With this interval division scheme we extend an approach presented in the pioneering 
paper by Lofstedt et al. (1993). In particular, we will treat the bubble collapse phase 
denoted by C in figure ^(6) in the interval tm < t < t* , the reexpansion interval (R) very 
close to the time of maximum compression {t « t*), the afterbounces (AB) for t* < t < 
and the bubble expansion in two stages for < t < (Ei) and < t < t,n (E2). 



Within each of these intervals, certain pressure terms in (2.1) are dominant, whereas 
others are negligible. Thus, simplified equations with analytical solutions can be derived, 
which enable us to characterize the complex bubble behaviour analytically and quanti- 
tatively. Our approximate formulas hold in the regime of strongly collapsing bubbles, 
i.e., for Rq > R*J'{Pa); in the weakly oscillating regime, the bubble dynamics becomes of 
course trivial. 

4.1. Rayleigh collapse (region C) 

We first take a closer look at the main collapse (after Rmax has been reached, interval C in 
figure ^j) . Figure ^ shows the behaviour of the most important terms in the RP equation 
(defined in table [l]) just prior to the main collapse. The abscissa displays the logarithm of 
the time interval before the collapse time t* which is identified by the condition R{t*) — 0, 
i.e., the bubble reaches its minimum radius at t* . The ordinate gives the logarithms of 
the absolute values of the various pressure contributions. As the whole time interval 
treated in this subsection only comprises « 0.1 /zs, and we want to discuss processes as 
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fast as Ins, we choose the polytropic exponent in (2.3) to be k = 5/3, the adiabatic 
value for argon. Note that the portfons of the graphs for \t* ~ t\ < IQ^'^T in figure ||, as 
weh as in figures H - |l^ below, represent time scales on or below the picosecond scale. 
As hydrodynamics breaks down here, this part of the computation will only be able to 
give a reasonable effective dynamics. We will take care not to draw physical conclusions 
from data in this range. 




Figure 8. Relevant pressure contributions according to the complete RP equation on a log- log 
scale before the instant of collapse t* . The bubble is driven at Pa = lAatm and its ambient 
radius is Rq — 4fj,m. 

In a large part of the collapse phase (figure H) the dynamical terms Pacc and pyei give 
the dominant contribution; they compensate each other, so that the dynamics is well 
described by the classical Rayleigh collapse 

RR+^R^ = 0. (4.1) 



This formula complements the quasistatic approximation (3.1) above. Equation 4.1 
plies a scaling law for R{t): 

2/5 



im- 



R{t) = Rn { *—^) . (4.2) 
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Here, the oscillation period T is used for non-dimensionalization of the time coordi- 
nate. The characteristic radius Rr can be estimated from an energy argument: at 
R — Rmax, the potential energy of the bubble is approximately Epot ^ AttPoR'^^^^/S, see 
e.g. Smereka, Birnir & Banerjee (1987). Converting this into kinetic energy of the fluid 



at i? = i?o, we get as an estimate for the bubble wall speed a.t R = Rq 



R 



R=Ro 



2Po 



1/2 



R„ 



Rq 



3/2 



Using the time derivative of (^j), we find 

2/5 

f25PoT^ 



5T 



Rr — Ro 



R\r= 



R=Ro 



2Rn 



1/5 



RfH^^U.Sfim- 



3/5 



(4.3) 



(4.4) 



With this Rji, (4.2) is compared to the numerical result of the RP ODE in figure y. 
Both slope and prefactor are reproduced excellently, despite the rather crude approxi- 



mations leading to (4.4). The only characteristic value for the Rayleigh collapse is Rmax, 
which depends on Pa and (although weakly) on Rq. Analytical expressions for these 



dependences will be given in Section 4.4 



We now examine the range of validity of (4.1); one could worry whether it is justified 
to neglect pgas and Psnd during collapse. For the solution ( |4.2| ), we have p^ei = —Pacc oc 
{t* - ^)"^^^ whereas (as long as R{tf > h^) Pgas oc {t* ~ t)-"^ and p,nd oc {t* - i)""/^ 
for K, — 5/3, i.e., the latter two pressure contributions grow stronger than the dynamical 
terms as t — > . This can also be observed in figure ||, but the absolute value of Pgas and 
Psnd is negligible compared to p^ei , Pace except for times very close to t* . We can compute 
the range of validity of (4.1) by equating pgas = Pvei and Psnd = Pvei, respectively, using 
(4.2), (4.4). It turns out that the sound pressure contribution is the first to violate ([4.l| ). 
This happens at tsnd with 

1/7 



Rn 



(1 



^5/7 



Ro 



Ro 

Rm.a: 



max 
18/7 



18/7 



ciT 

« l.Ox 10-'*(l + as)^/^ 
which agrees with the numerical result e.g. in figure |^ (where i?o = 4 pim and i?„ 



Ro\ 
lim J 



(4.5) 



47 fim). For the approximation ( }4.5D , R'^{tsnd) ^ was assumed; is the surface 



tension parameter introduced in (2.6). The collapse behaviour changes due to Psnd shortly 
before another assumption for (4.1) breaks down: obviously, R{t) cannot be smaller than 
the van der Waals hard core h. Equating R{t) — h using (4.2) with h — _Ro/8.86, we 
obtain the "hard core time" 



{t* — tvdw)/T 



\25Po 



1/2 



1 /l5/2 



tr 



,3/2 



5.5 X 10" 



^0 

Rmau 



3/2 



Rq 



^ ■ (4-6) 



At i ~ tydwi the van der Waals hard core cuts off the scaling behaviour abruptly. However, 
for typical values of i?o ~ 4 fim, the bubble collapses like an empty cavity for a time 
interval from {t* - t) ^ 1 ps down to {t* - t) ~ 100 ps {t* - t ^ 0.03T . . . 3 x lO^^T). 
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(t*-t)/T 



Figure 9. Collapse dynamics of R{t) for the same parameters as in figure g. This bubble reaches 
a maximum radius Rmax ~ 47 fim before collapsing. The theoretically expected Rayleigh scaling 
law R{t) = RR{t* — t)^^^ (dashed) is followed accurately. Also indicated are the limiting times 
tvdw,tsnd for the vahdity of the scaling law. 



4.2. Turnaround and delayed reexpansion (region R) 

As the gas is compressed to the hard core radius, the collapse is halted abruptly. Lofstedt 
et al. (1993) have shown that - in the Hamiltonian limit neglecting PsunPvis , Psnd and the 
temporal variation of p^xt ~ the turnaround time interval of the bubble is approximately 

1/2 



Till 



R{t*] 
R{t*] 



1/2 



Po(l 



h 
Rq 



3k 



Rn 



(4.7) 



This equation also follows from approximating the RP equation ( 2T ) by keeping only the 
dominant terms in the immediate vicinity of the collapse, i.e., Pacc and pgas (cf. figures ^ 
and |l^): 



piRR = Po (1 + as) 



\R^{t) - 



(4.8) 



(4.8) is a good description of bubble dynamics for a time interval around the collapse of 

length ~ Tturn- 

Figure O shows the pressure contributions after t* . From {t — t*)/T ~ 10^^ to 
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Figure 10. Important pressure contributions according to the RP equation on a log- log scale 
after the instant of collapse. The same parameters as in figure H apply. 



{t - t*)/T - lO""* (i.e., from {t - t*) - 30ps to ~ 3ns) the dominant terms in 
are pgas and Psnd, which compensate each other. This means that the energy stored 
in the compressed gas is released almost exclusively through emission of sound waves 
(cf. Church 1989) - it is not converted back to kinetic energy of the liquid surrounding 
the bubble. The corresponding dynamics shows a relatively low expansion velocity and 
small acceleration, keeping a very small bubble radius for a few ns (figure pi] ) . This time 
interval of delayed reexpansion (denoted by R in figure ||6) is described by 

R{t) d 



0=Pgas(R{t)) 



ci dt 



PgasiRit)) 



with Pgas given by (2.3). This ODE has an analytical solution 



6k, 
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(4.10) 



For {R{t) - Rmin) < R 

simplified to yield 



R{t) K R. 



i.e., just after the collapse, this implicit equation can be 
c, R^ . — 



3k Ri 



(4.11) 



22 S. Hilgenfeldt, M. P. Brenner, S. Grossmann, and D. Lohse 

This linear expansion law holds for a longer time interval if Rmin is larger, i.e., for smaller 
Pa ■ Its validity is demonstrated in figure |ri|. Note that although the turnaround time 
Tturn becomes smaller for decreasing Rmin — h, the velocity of the bubble immediately 
after collapse is actually smaller because of the larger energy losses through acoustic 
radiation. The strongly asymmetric shape of R{t) around t* has also been observed 



in experimental measurements of bubble dynamics, e.g. by Barber & Putterman (1992) 



Tian et al. (1996), [Weninger, Barber fc Putterman (1997)|, and|Matula et al. (1997)|. 
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Figure 1 1 . Bubble dynamics for 7?o = 5 /im and two different driving pressure amplitudes 
p = 1.3 and p = 1.5 in the vi cinity of collapse. The dashed lines give the linear reexpansion 
approximation from equation ( 4.11 ). 



After the delayed reexpansion phase, the bubble wall gains speed and enters another 
short time interval around {t — t*) ~ 10~^T well described by Rayleigh's equation (4.1) 
with R{t) cx (t — t*)^/^ as the bubble expands. At (t — t*) « 10~^r it enters the phase 
of subsequent afterbounces. 



4.3. Afterbounces: a parametric resonance (region AB) 

The discussion of the afterbounce interval (AB in figure ^) is intimately connected to 
the explanation of the wiggly structure of various dynamically computed terms, like 
the expansion ratio (figure ^6) or the diffusive equilibrium lines in Fig ^ Obviously, 
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as the RP equation (2.1) describes a driven oscillator, the maxima in the expansion 
ratio represent parameter values of resonant driving. Figure |lj clarifies the character 
of these resonances. It shows two time series of the bubble radius R{t) at values of Rq 
corresponding to a relative maximum and a relative minimum of Rmax/Ro, respectively. 
A large or small expansion ratio results from the phase of the afterbounces at the time 
when Pf^j-t becomes negative, i.e., when the external forces start the rapid expansion: for 
the bubble with the large expansion ratio, the last afterbounce "fits" into the expansion, 
which is enhanced. For the other bubble, growth is inhibited as the last afterbounce 
collapse interferes with the expanding external force. 
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Figure 12. Two time series R{t) for Pa = 1.2 atm and Tio = 7.05 /im (dashed), Ro = 7.30 /im 
(solid). The first ambient radius corresponds to a resonance minimum in Rmax/Ro, the second 
to a maximum. Note that the afterbounces at the beginning of the main expansion are just one 
half-cycle out of phase. 



The afterbounce oscillations show relatively small amplitude, and it is therefore pos- 
sible to linearize the RP equation in this region of the driving cycle. Moreover, sound 
radiation and viscosity contributions are negligible. For simplicity, we also neglect Psur 
for the moment. In order to separate the time scale 1/uj of the driving from the much 
shorter time scale of the afterbounces, which is ~ 1/we, we use the ansatz (cf. e.g. 
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Hinch 1991) 

R{t) = Ro{T)il + y{t)) (4.12) 

with small y(t) and a slowly varying function Ro{t) which is to be determined; r = 
et, e = uj/uje ^ 1- To leading orders in y and e, equation (2.1) is transformed into 



- — (1 — p COS LUeT) 



(4.13) 



Requiring the slowly varying (secular) term on the right-hand side to vanish, we have to 
choose 

Ro{r) ^ Ro/ {I - PCOs{lu,t))^^'' ^ Ro/ {1 ~ pcos{Lot)f . (4.14) 

(4.15) 



With this definition, ( 4.13 ) results in a Hill equation: 

y + Wg(l - pcosutf/^y = 0. 



Because of the separation of time scales S> w this equation represents a harmonic 
oscillator with slowly varying eigenfrequency, i.e., the afterbounce frequency uJab = 
We(l — pcosujt)^/^ . For this system, E{uJab)/i^ab (with E = (2/^)aj^j/2 being the os- 
cillator energy) is an adiabatic invariant (see Hinch 1991), i.e.. 



{y^){l -pcosujtf/^ = const. , 



(4.16) 



where the mean ( • ) is an average over the fast time scale. Note that in the time interval 
7r/2 < Ljt < 37r/2 of aftcrbounces (1 —pcosujt) > 0. It follows that the amplitude of the 
afterbounces changes as Roy oc (1 — pcosuit)^^^'^. 

The resonance structure of ( 4.15 ) still cannot be evaluated analytically. Yet the para- 
metric driving of (4.15) has a very similar shape to the cosine driving of a Mathieu 
equation. We can therefore further approximate ( 4.15| ) by choosing suitable constants 
Qi,Q2, where we require 



Qi ^ Q2 cos{ujt) = (1 — p cos{u!t)) 



5/3 



for ujt — —,TT . 

2 ■ 



(4.17) 



i.e., Qi = 1, Q2 = (1+p)^^'^— 1. The errors in this approximation are only a few percent in 
the time interval ^ [7r/2, 37r/2] of afterbounces we focus on. As an analytically accessible 
approximation to the afterbounce dynamics of (2.1) we have thus the Mathieu equation 



y"+44 (I-Ul+P)'/' 



1 



COS 2x\y = 



(4.18) 



with dimensionless time x = u)t/2] the primes denote derivatives with respect to x. 

The contribution of surface tension may be included if <C 1 to yield a refinement of 
(HH): 

y" + ^4^(^{l+'^-a)j- (l+p)5/3-l + ^a, (l + 2p-(l+p)5/3) cos2i^y = 0, 

(4.19) 

ojs does not cover 



with the factor as from (2.6). Note that a simple substitution Ue 
all first-order effects of 
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For certain parameter combinations, equation (4.19) shows parametrically stable or 
unstable solutions. Because We/w 3> 1, the best analytical approximation to these char- 
acteristic values is given by the asymptotic series (Abramowitz & Stegun 1972) 

b^.V-s (4.20) 

with b = 4^[(l-2a,/3)(l+p)5/3+4a,/3-(l+p)], s = 8^[(l-2a,/3)(l+p)5/3+4a,/3- 
p — 1 + 2as/3], and v — 2kM ± 1, where the sign distinguishes even from odd Mathieu 
solutions. kM is the order of the Mathieu resonance, corresponding to the number of 
afterbounces in the RP equation (see below) . We take here only the leading term on the 
right-hand side of ( 4.20| ) and treat the case kM ^ 1, so that v « 2kM] moreover, we 
only keep terms up to first order in a^. This yields ambient radii Rq^''''^ for which the 
oscillation shows maximum stability against parametric excitation: 

r,ik.,)^f 3Po q'/' J_ 2a / g5/3_2g+l 2 - gV^ 

« 74.0 Mm -^^=-^ + 0.487 Mm - ^_^i±i + 2^^^ ) (4.21) 



Vg5/3 _ 1 kM ' ■ ^ y q5/3_i q2/3 

Here we have abbreviated q = (p + 1). Note that the correction term due to surface 
tension does not depend on /cm- 

Although the behaviour of the RP oscillator is fairly well described by Mathieu os- 
cillations in the afterbounce phase (see e.g. figure |l^), it is of course entirely different 
during the expansion interval of the cycle. Therefore, some information about the overall 
shape of the oscillation must enter into our analysis. Especially, Mathieu solutions can 
be T- or 2T-periodic. RP dynamics in the SL regime, however, only allows for T-periodic 
solutions, as the 2T-periodic Mathieu solutions would require large negative values for y. 
Therefore, every second resonance of ( 4.19D must be dropped, i.e., the resonance of order 



kM of ( 4.19 ) corresponds to resonance number k = kM/'2 of (|2.l| ), so that the k*^ reso- 
nance radius ^Rq*^'* of (4.19) for a dynamics with k afterbounces is obtained by replacing 
kM by 2k. 

We must also provide additional information about the length of the afterbounce inter- 
val. Figure |l^ shows that this length, which is almost independent of p for our Mathieu 
approximation, is significantly reduced for increasing p in the case of the RP equation, 
as the expansion interval lasts longer. This is a property of the nonlinear part of the 
RP cycle (cf. the next subsection), which cannot be modelled within the Mathieu ap- 
proximation. In the Appendix it is shown that the length of the afterbounce phase (and 
therefore the resonance number k) has to be rescaled according to fc ^ C{p) x k with C{p) 
approximately given by the expansion C{p) « 0.688 — 0.548(p — 7r/2) + 0.418(p — 7r/2)^ 

(cf. (ii)). 

Figure |lj presents a comparison of the computed resonance radii of order (number of 
afterbounces) k from numerical solutions of ( ^.l[ ) - both for relative maxima and minima 
of Rmax{Ro) - and from (4.1£) for driving pressure amplitudes p = 1.2 and 1.5, corrected 
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Figure 13. Comparison of time series from the numerical solution of the RP equation ( |2.l[ ) 
(upper) and a T-periodic solution of the Mathieu approximation ( [4.19| ) (lower) for the same 
parameters p = 1.5, -Ro ~ 6.0 /im. The onset of afterbounces is delayed in the full RP dynamics. 
The correction factor (B2) has to be employed. As the Mathieu solution is divergent (grows 
exponentially from cycle to cycle), the lower curve was normalized once per cycle. 



with C'{p). The resonance locations are in good agreement for both pressures, considering 
the multitude of approximations they were calculated with. 

(k) 

In the stability maxima marked by Rq the bubble is less excitable by the driving 
than bubbles with neighboring Rq. Therefore the expansion ratio has a local minimum 
and the average pressure {pgas)^ oc l/Rf^^^^ experiences a local maximum. The existence 
of such wiggles in Rmax/Ro (and therefore in (pgas)^ leads to the possibility of multiple 
equilibria for given experimental parameters Pa and poc i the ambient radius can adjust 
itself diffusively to different stable equilibrium values, depending on initial conditions 
and/or perturbations. For an analysis of the physical consequences of multiple equilibria 
we refer the reader to the work of Brenner et al. (1996a), Hilgenfeldt et al. (1996) and 



Crum fc Cordry (1994) 
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Figure 14. Ambient radii R^^'^ as a function of resonance order k = fcM/2 from the full solution 
of ( |2.l[ ) (circles) and the analytical Mathieu approximation (4.21) (crosses); Pa = 1.2 aim (left) 
and 1.5atm (right). Note that the circles represent relative maxima and minima of Rmax, 
because the Mathieu equation admits twice as many resonance values as the RP equation. The 
Mathieu solutions were rescaled by the factor C{p) of equation (B2). 



4.4. Bubble Expansion (regions Ei, E2) 

We wish to be more quantitative about the properties of the expansion phase now. 
Despite the small portion of parameter space for SL bubbles, there are different types 
of expansion behaviour to be identified depending on p and Rq. For p ^ 1 and large 
Rq > 10 /im, the gas pressure plays an important role and balances the dynamical pres- 
sure, which is dominated by p^jg; for most of the cycle, so that a first approximation to 
the dynamics is: 

3 



Pgas{R,t) (4.22) 

With pgas{R,t) ~ PqRq/R^ for large Rq, this equation yields a solution for R(t): 



R 



5/2 



2^ 

3/3/ 



1/2 



2/5 



(4.23) 
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with the starting time of expansion t_ and starting radius R- — R(t-). For radii 
R 3> R-, ( 4.23 ) reduces to a Rayleigh-type expansion law, which gives a scahng relation 
for the maximum radius: 

Rmacc (X RI^^ (4.24) 

if we assume that the length of the expansion interval is independent of Rq, which is a 
good assumption except for small Rq below the Blake threshold {iA). The law ( 4.24 ) 
can numerically be confirmed for large Rq, see Hilgenfeldt et al. (1996). Together with 
( ^■llD this yields (pgas)4 oc rI'^ . 

For higher driving pressure amplitudes or smaller Rq, i.e., in most of the SL parameter 
region, the approximation ( 4.22 ) is too crude. Instead, one has to take into account the 
balance of the dynamical pressures PaccPvei and the external pressure p^xt 

pi(^RR+^R^^ ^Po{p cos Lut-1), (4.25) 

as can be seen e.g. in figure 




Figure 15. Upper: important pressure contributions for Pa — lAatm, Rq — 4.0 nm during 
the expansion phase. The thick dashed line gives the external pressure pext = Po — PaCosut. 
In the time interval Ei, it is primarily balanced by Pvei (solid), in E2 by pacc (dashed). Lower: 
bubble radius expansion dynamics for the same parameters. 
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In the work of Lofstedt et al. (1993), the left-hand side of (4.25) has been approximated 
using a power series for R{t). This leads to a bubble expansion which is linear in time, 
with a velocity proportional to y/p — 1. The first nonlinear corrections are of fourth order 
in t. 




Figure 16. Comparison of RP dynamics (solid) and theoretical approximation (dashed) of 
bubble expansion dynamics for Pa = 1.5 atm, Ro — 5.0 nm. The theoretical solution becomes 
complex left and right of its zeros. Also shown is the suggested theoretical solution (thin dashed 
line) from equation (29) of Lofstedt et al. (1993), where R{t = 0) was matched to the value of 
the RP dynamics. 



However, an expansion like this turns out to be a series which is not well-controlled. 
Especially, no reliable values for the time tmax and value Rmax of the maximal radius can 
be derived (cf. figure ^6|). We therefore follow a different ansatz: consider the dynamical 
terms Pacc and Pvei for a typical bubble expansion (figure ^|). During the first, almost 
linear part \pvei\ S> |pacc|, whereas when the maximum is approached and the bubble 
decelerates, \pacc\ 3> Ipueil- This suggests a division of the expansion interval into two 
parts (denoted by Ei and E2 in figures §(&),|l^, and p^. Observing that the combination 
RR + R? is just the second derivative of i?^/2, (4.25) can be approximated with good 
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accuracy by 

d^ 

dx"^ 
d^ 
dx^ 



R'' = -Ri^^ip cos X 



4. 
9^ 
2 

— J 

3 



R^ = -RL.ip cos X 



for 



for 



< X < Xa 



(4.26a) 



(4.26b) 



Here we have introduced the dimcnsionless time x = cot and the hnear bubble resonance 
radius Rres from (|2J); Xm is given by i?(x,„) = Rmai 
right-hand side make sure that the dominant terms p^ei in (4.26a) and Pa 



The rational prefactors on the 
in (4.26&) are 



correctly represented, while the other terms gives contributions of the indicated order in 
each case. The starting time a;_ = a;_(p) and the transition time a;+ — x+{p) between 
both solutions are given by the zeros of Pe^^t, i.e., x^{p) = —x^{p) = arccosl/p. Equa- 
tions ([4.26|a,6) can be integrated analytically requiring continuity and differentiability 
at x+ for the overall solution. To complete the problem, initial conditions at X- have 
to be imposed: we set i?_ = i?(a;_) ~ (^Rq with a parameter C ~ 1 whose value is not 
crucial for the shape of the solution. An estimate of C, can be computed from algebraic 
equations, but not in an explicit form. C. lies between 1.2 and 2.0 for typical Rq of SL 
bubbles; for simplicity, we choose C = 1-6 in all calculations. For the initial velocity, 
we observe that X- marks the transition from the afterbounce regime, where the bubble 
essentially oscillates with its cigenfrequency, to the expansion regime, where the govern- 
ing time scale is the driving period T. Therefore, we set R'_ = (dR/dx){x^) = 
corresponding to i?_ = luR^ in dimensional terms. 

Figure shows that the shape of the expansion as well as time and value of the 



maximum are reproduced satisfactorily. From the solutions of (4.263) in Ei and (4.266) 
in E2 one obtains a system of equations for Rmax and x^'- 

RLax=R'- (l + 2(a;„ + a;+)) 



-R^ 



I — p cos x„ 



psmxr, 



+ i(psinx+ - x+){x,, 



I, , 3i?i 

+ -(psma;+ - x+) + = 0. 



3X4 



,(4.27) 



(4.28) 



Note that ( 4.27 ), ( 4.28| ) give the position and height of the radius maximum without any 
freely adjustable parameters. The inset of figure |l^ shows the maximum radii obtained 
with these formulas for i?o = 5 /im and 9 /im and p = 1 — 5 together with results from 
a complete RP computation. Apart from the resonance wiggles (cf. Section ^^), the 
curves are very well reproduced both within the p regime for SL bubbles and for the 
much higher pressure amplitudes of cavitation field experiments. 

Equation ( |4.2^ ) is transcendental, and Rmax and Xm do not have simple analytical 
representations. One can, however, derive simple expressions in several limiting cases: if 
p 3> 1 , we obtain after a lengthy calculation 

« y^dRrUFp-Gf/^ 



R-n 



(4.29) 



with constants F = l-t-57r/6 = 3.618 . . . , G = 197rV24- 1 = 6.813 .... This formula is a 
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Figure 17. Comparison of direct numerical results of Rmax (solid) for _Ro ~ 5 (lower) and 
Ro = 9 j-im (upper) in the p ran ge relevant for SL with the approximate law (4.31) (dot-dashed) 
and the simp ler fo r mula ( 4.3S ) (long dashed). The inset shows that the full theory (dashed) 
according to ( 4.27 ), ( 4. 26 ) for Rmax{p) for the same two Ro gives excellent agreement with the 
direct RP computation (solid) for both the SL range and the cavitation field regime p ~ 2 — 5. 
Even the approximation (4.31) yields good results up to p ~ 3 (dot-dashed). 



good approximation only if p > 5. For p > 1, i.e., the case of interest for sonoluminescent 
bubbles, we can expand Xm around x = 7r/2. Moreover, as Rres 3> i?o for SL driving 



frequencies, we can also neglect the last term on the right-hand side of (4.28). To leading 
order in {xm — 7r/2), ( 4.28 ) then becomes 



i(psina;^ 



x+) ; 



(4.30) 



remember that x+ — x+{p) — arccos(l/p). For p e [1.0,1.5], the second term of the 
right-hand side of this equation is never greater than 0.185p in absolute magnitude, so 



that good approximation. Inserting into ( |4.27| ) gives 

Rliax = fip)Rl + 9{p)Rles with 
/(p)-C'(l + 2(p + x+)) , 



9{p) = 



1 -pcosp- ^(jp- - x\) + ^{ps\nx+ - x+){p + 2,x+) 



(4.31a) 
(4.31b) 

(4.31c) 
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The second term in (4.31a) is much greater than the first; thus, it is not primarily 
Ro which determines Rmax, but the resonance radius Rres- With Rres oc 1/w, we see 
that the expansion ratio is (at constant p and Rq) roughly inversely proportional to the 
driving frequency, i.e., upscaling of SL collapses can be achieved by lowering uj, which 
was also seen in experiment by R. E. Apfel (private communication, 1996). In the same 
way, a higher ambient pressure Pq (while keeping p — Pa/ Pq constant) will lead to 
higher expansion ratios because of the dependence of Rres on Pq (see (^)). A further 
simplification of (4.31) can be obtained from a stringent expansion in [p— 7r/2) and Ro, 
which yields to leading order the simple result 



R-n 



= 67.2 + 0.112 



^) +99.5(p 
fim 



n/2) + 0{{p-n/2f 



(4.32) 



Figure shows the very good agreement of (4.31) and ( 4.32 ) with full RP dynamics 
for several i?o over the whole range of pressures in SL experiments. The approximation 
breaks down only at p ~ 3, where Xm ~ p is no longer valid, see inset of figure ^ The 
expansion ratio is also accurately reproduced for moderate or large Rq by this formula, as 
seen in figure |l^. The deviations for small Rq are due to neglecting Psur, which becomes 
the dominant influence as Rq approaches R^. 

One would therefore like to include the effects of surface tension into (4.31). We make 
the following ansatz: instead of ( 4.26| a, 5), we write 



R' 



dx"^ 
da;2 3 



p cos X 



p COS X 



1 



K{p) 
K{p) 



for 



for 



< X < X4 



(4.33a) 
(4.33b) 



This models the influence of Psur by an average pressure contribution of PoOs/ K{p), 
where K is taken to be independent of Rq. Expanding Xm again around 7r/2, we get 



Pfnax = f{p)Rl 



gip) 



3K{p) \2 



1 

r^px 



Rr- 



with /(p), g{p) from (4.316, c). With this expression, (i?„ 
maximum at 



Rl{p) 



3a W 



\px+ 



(4.34) 

/Ro){Ro) shows a global 
(4.35) 



PoKip) g{p) 

For large enough p, we can equate ( [4.35|) and (3^), because R\^ and Rq are very close 
then. This gives an estimate for K{p): 

\{p^ + x\) + \px+ 



K{p)^-^Vi{P'l) 



gip) 



(4.36) 



In the regime of SL driving pressures (p = 1.2 . . . 1.5) K{p) depends only weakly on p; 
its value lies between 7.5 and 9.4. The ansatz ( 4.33| n, b) proves very successful for the 
description of {Rmax / Ro) (Ra) over the whole range of relevant Rq, as can be seen from 
figure |8|. 

From (4.34), we obtain expansion ratios and, using (2.9) and (2.11) for given gas 
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Figure 18. Direct numerical computation (solid) and theory (dot-dashed) according to (4.31) 
for th e exp ansion ratio {Rmax / Ro){Ro) at p — 1.2 (lower) and 1.5 (upper). The dashed curve 
from (4.34) takes surface tension into account and is able to reproduce the numerical graphs for 
almost all i?o. 



concentration Poo/ Po, an approximation for the location Ro{p) of diffusive equilibria can 
be computed. Figure [l^ shows that both the stable and the unstable branches of the 
equilibrium curves (taken from figure |l|) are reproduced satisfactorily for both high and 
low gas concentrations. 

When, starting on the stable branch, p is lowered, Rmax/Ro becomes smaller and, by 
( p.ll ), {pgas)i becomes larger. The corresponding equilibrium ambient radii i?o shrink. 
Eventually, the minimum of {pgas) 4^/ becomes larger than Poo/Po (see figure^) and no 
Rq can fulfill the equilibrium condition. This situation corresponds to the turning point 
of the diffusive equilibria in figure |l| and figure |l^. 



5. Role of surface tension and liquid viscosity 

The previous sections have provided a detailed analysis of the dynamics of SL bubbles. 
How will these results change if we introduce a different fluid with different surface tension 
a and/or fluid viscosity r/i? 

Surface tension is the crucial parameter for the location of the Blake threshold in 
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Figure 19. Diffusive equilibrium lines from figure |l] (solid) and their approximation using 
(4.34) (dashed) near the parameter domain of sonoluminescing bubbles. 



parameter space (cf. also Lofstedt et al. 1995 or Akhatov et al. 1997). The transition 
from weakly oscillating to strongly collapsing bubbles and therefore the boundary of the 
SL region determined by (1.1) is entirely controlled by a. If we had ct — > 0, bubbles 



with any i?o would be strongly collapsing, i.e., liquids with small surface tension should 
allow for violent collapses at smaller Pa- On the other hand, in liquids with high a larger 
Pa and i?o are required for collapses. Although a larger a has a stabilizing effect on 
the bubble surface, 



^ and ( 4.35 ) show that the \Mg\ — 1 line overtakes the shape 
instability threshold in Pa-Ro parameter space (cf. Hilgenfeldt et al. 1996), so that no 
stable SL should be possible if a is e.g. five times higher than in water. This is easily 
confirmed by the numerical solution of the RP equation, see figure pQ(a) . 



At first sight it seems that fluid viscosity could have been neglected in (2.1) right from 
the start. Apart from a slight damping effect during the afterbounce phase, the influence 
of rji for water on bubble dynamics is hardly noticeable, even a tenfold increase of r/; only 
reduces the maximum radius by w 10% (figure ^l|a,6). We can, however, estimate by how 
much the viscosity would have to be enhanced to have a significant effect: the damping 
of a high viscosity liquid should ultimately prevent the bubble from collapsing violently 



and therefore it will never fulfill the energy transfer condition (1.1). As the collapse is 
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Figure 20. Parametric stability thresholds (solid), calculated as elaborated in Hilgenfeldt 
et al. (1996), and \Mg \ — 1 lines (dashed) for bubbles of different surface tension (a) and viscosity 
(fe). All other parameters were fixed at the values for water. Increasing a makes the region of 
potential SL bubbles above the \Mg\ = 1 and below the stability threshold vanish, whereas 
higher rji enlarges this area. 



in fact the first afterbounce minimum, an estimate for this critical rjf can be obtained if 
we demand the afterbounces to be overdamped. The viscosity r]i introduces a damping 
term in the linearization of the RP equation (2.4). It is easy to see that overdamped 
motion requires 



With the definition of LUe 



2^F 
it follows 



> 



1/2 



(5.1) 



(5.2) 



With a typical value Rq = Afim and keeping a, pi at the values for water, we obtain 
that rjf > 4:0r]„ater- This is confirmed by direct computation of ( p.l[ ) using r]f and strong 
driving pressure Pa — lAatm, see figure ^(c)-(e). Viscosities in this range can be 
easily achieved in mixtures of water and glycerine. For moderate glycerine percentage, 
the viscosity is not very different from pure water, but for high glycerine contents it 
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rises dramatically. The required factor of 40 is (at 10° C) reached for « 70% glycerine 
(weight percentage). Above this percentage, it would be extremely difficult to obtain 
collapses strong enough to ensure energy transfer and the ignition of SL. Moreover, 
chemical dissociation reactions in air cannot take place, which seem to be necessary for 
SL stability using air at moderate degassing levels (Lohse et al. 1997). This may be 
the reason why Gaitan (1990) was not able to observe stable bubbles above a glycerine 
percentage of « 60%. The actual threshold for SL should occur at slightly smaller rji than 
predicted by (5^), because even if the collapse minimum is not completely damped out, 
the collapse is already considerably weakened (figure pl|c). Also, the threshold should 
be higher because of the additional damping effect of thermal dissipation (see Vuong & 
Szeri 1996, Yasui 1995) which is not included in our approach. 
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Figure 21. Bubble dynamics R{t) for diflferent viscosities r]i/r)^ater ~ !(«), 10(6), 35(c), 40(d), 
45(e). Ro = 4/^m and Pa — lAatm in all cases. The transition from collapse / afterbounce 
behaviour to aperiodic dynamics occurs near rfi/rj^ater =40, as predicted from (5.2). 



Even for smaller rji ~ rj^ater , fluid viscosity is an important contribution to the damping 
of bubble surface oscillations, as was shown in Hilgenfeldt et al. (1996). Therefore, a 
moderate increase in rji does not lead to significant changes in the R{t) dynamics itself 
(and therefore in the \Mg\ = 1 curve), but it helps to stabihze bubbles at larger Rq. 
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This change affects only the parametric and afterbounce instabihties (see Hilgenfeldt 
et al. 1996), which can show accumulation effects over several driving periods, but not the 
Rayleigh- Taylor instability, which is directly dependent on the acceleration of the bubble 
wall and cannot change significantly when the R{t) dynamics does not. Therefore, the 
Rayleigh- Taylor process still cuts off the bubble stability region (cf. figure 6 of Hilgenfeldt 
et al. 1996) at Pa ~ lA5atm almost independent of viscosity, whereas at smaller Pa, 
much larger bubbles can be stable if rji is enhanced. 

In figure ^(6), we show surface stability curves for different fluid viscosities. For high 
rji, the region of stably sonoluminescing bubbles between the almost unaffected \Mg\ — 1 
curve and the increased stability threshold is considerably enlarged. This would probably 
correspond to a substantial upscaling of SL intensity. Gaitan's (1990) experimental 
observation that a moderate percentage of glycerine helps to establish stable SL bubbles 
supports this conjecture (see also the experimental results of Gaitan et al. 1996 for fluids 
of different viscosity and surface tension) . Combining our results for a and rji , we conclude 
that the ideal fluid for violently collapsing, but surface stable bubbles should have small 
surface tension and high viscosity. 

6. Conclusions 

The analysis of the RP equation presented here has explained quantitatively quite 
a lot of features seen in numerical computations of RP dynamics. One cycle of os- 
cillation of this highly nonlinear system can be completely divided into subsections in 
which its behaviour can be accurately approximated by analytically integrable equations. 
While being in the spirit of Lofstedt et al. 's (1993) previous analysis, the present work 
presents more complete and more detailed results. We emphasize that there are no freely 
adjustable parameters in our approach. 

We made use of these approximations to calculate analytical laws for the bubble's 
collapse, afterbounce behaviour, expansion dynamics, maximum radius, and expansion 
ratio. With these results we could clarify parameter dependences of numerically calcu- 
lated curves of diffusive equilibria in Pa-Ro parameter space as those in figure ^ A 
summary of relevant analytical relations and predictions for experimental verification 
was already given in the Introduction. 

An approximation of the RP equation by a Mathieu equation has explained the wig- 
gly resonance structure characteristic for many quantities derived from RP dynamics. 
The concept of a quasistatic Blake threshold between regimes of weakly oscillating and 
strongly collapsing bubbles was able to shed light on the existence of stable diffusive 
equilibria in the SL regime for high driving pressures. The change of sign in the slope 
of ipgas) ^ (Ro) is a generic feature of RP dynamics, resulting from the dominance of 
surface tension pressure at small Rq. This allows the bubble to reach a stable diffusive 
equilibrium. 

In all approximations of the RP equation, the fluid viscosity term for water or similar 
liquids could be neglected without causing large errors. Both numerical computations and 
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analytical estimates of the magnitude of pvis show that rji has to be as high as ~ 40 rj^ater 
to become a dominant contribution to bubble dynamics. Viscosity does, however, have 
a strong influence on the dynamics of surface oscillations; parametric instabilities are 
weakened for larger rji. 

Surface tension is the underlying cause for the change from unstable to stable diffusive 
equilibria, which stabilizes small bubbles to an extent that they can only show weak 
oscillations. For fluids with low tr, collapses of bubbles of a given size are more violent. 
This is especially interesting if this effect is combined with higher fluid viscosity to 
establish bubbles which show a similarly violent collapse dynamics as bubbles in water 
while maintaining larger radii. 

Other possibilities for an upscaling of the collapse intensity are the use of lower driving 
frequencies uj or of larger ambient pressures Pq at the same Pa/Po- These predictions offer 
a useful guideline to experimenters in search of upscaled single bubble sonoluminescence. 

This work has been supported by the DFG uder grant SFB 185-D8 and by the joint 
DAAD/NSF Program for International Scientific Exchange. 



Appendix A. Modifications of the RP ODE 



As stated in the Introduction, a lot of variations to the RP equation (2.1) are known 



from literature, see Lastman & Wentzell (1981, 1982) for an overview. We mention here 
the form derived by Flynn (1975a, 1975b) 



Pi 



RR 




PgasiR.t) - P{t) - Po - 



Pgas{R,t) ~ '^r^l—- — 



2a 



(Al) 



which contains correction terms of higher order in Rjci. It also includes time derivatives 
of the surface tension and viscosity terms. 

Gilmore's equation (see e.g. Hickling 1963) differs from the other RP variations in that 
its key variable is not pressure, but the enthalpy H of the gas at the bubble wall: 



Pi 



1 



RR 





R d 
G dt 



U. (A 2) 



Here, the sound speed C; is not a constant, but depends on H. The exact form of 
this dependence has to be specified by an equation of state for water, e.g. of modified 
Tait form (Prosperetti & Lezzi 1986, Cramer 1980). Gilmore's equation was shown in 



Hickling (1963) to be an accurate description of a collapsing cavity up to Mach numbers 



|M/| = \R\ICi as high as 5. Its validity for the present problem of collapsing gas bubbles 
is, however, not well established (Prosperetti 1984). 

Figure ||(a) compares the bubble radius dynamics R(t) computed from Eqs. (A 1) and 
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(A 2) to the solution of (2.1). On this scale, the curves are almost indistinguishable. 
Only a blow-up of the region around the radius minimum reveals deviations. 



As all equations , ( A 1 ) , ( A 2 ) have a common limit for small Mi , significant devi- 
ations are only to be expected during collapses, when the bubble wall velocity becomes 
of order of the sound speed. Figure ^(c) shows M; for the different R(t) dynamics. Ob- 
viously, large differences occur only around the main collapse, and they are not impor- 
tant for the overall dynamics of the bubble, the collapse time interval being exceedingly 
small. Thus, the RP equation provides a relatively simple and very accurate descrip- 
tion of bubble wall motion, which is confirmed by comparison with recent experimental 
measurements of R{t) in Tian et al. (1996). Note also that our criterion (1.1) for energy 
transfer uses Mg, the Mach number with respect to Cg, which (for argon) is almost 5 
times smaller than c; (in water). Therefore, the energy transfer threshold can be well 
computed within the RP-SL approach. 



Appendix B. Length of afterbounce interval 

The Mathieu model equation ( [1.19| ) can only be expected to give an accurate descrip- 
tion of RP dynamics in the time interval [7r/2, 37r/2], for which it was matched to the 
Hill equation ( 4.15 ) via ( 4.17 ). Therefore, we have to compare the afterbounce interval 
of RP dynamics to this constant interval of length tt. Figure |l^ shows that for large 
driving p the Mathieu dynamics gives a good approximation for the end point of the 
afterbounce interval (starting point of expansion phase), but fails to model the onset 
time of afterbounces, i.e., the dimensionless collapse time x* . Thus, the length of the 
afterbounce interval is smaller than tt by a factor of C{p) — (37r/2 — x*{p))/-K. 

The collapse time x* is relatively close to the maximum time x^- It is therefore 
convenient to compute it from an expansion of (4.266) in powers of (x— 7r/2) and (p— 7r/2). 
This yields 

X* « 2.55 + 1.72(p - 7r/2) - 1.31(p - n/if . (B 1) 

The coefficients can be computed analytically, but are of very complicated form. This 
expression corresponds to a correction factor for the length of the afterbounce interval 
and, equivalently, for the resonance order fc, of 

C « 0.688 - 0.548(p - 7r/2) + 0.418(p - . (B 2) 

This is the correction introduced in Section [4.3| which is important for a satisfactory 
description of the resonances of (2.1) by (4.19), see figure |l3. 
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